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ABSTRACT 


Thi3  paper  presents  a  method  of  obtaining  one  type  of  asymptotic 
solution  for  the  linear  hyperbolic  second-order  partial  differential  equa- 
tion in  n-1  independent  variables  and  one  time  variable  with  coefficients 
independent  of  time.  Appropriate  initial  and  boundary  conditions  are  speci- 
fied.  In  particular  it  is  supposed  that  a  time-harmonic  forcing  term  is 
present  either  in  the  non-homogeneous  term  of  the  differential  equation  or 
in  a  boundary  condition.  Then  the  "steady  state"  part  of  the  exact  solu- 
tion has  the  same  time  behavior.  The  spatial  behavior  of  the  "steady  state" 
part  is  then  expressed  asymptotically  as  a  series  in  positive  integral 
powers  of  —  ,  where  to  is  the  circular  frequency  of  the  source.  The  suc- 
cessive coefficients  of  this  series  are  respectively  discontinuities  of  a 
particular  solution  of  the  differential  equation  and  of  successive  time- 
derivatives  of  this  solution. 

By  a  study  of  discontinuities  of  solution: of  the  differential  equa- 
tion it  is  shown  that  the  discontinuities  needed  for  the  successive  coeffici- 
ents of  the  expansion  can  be  obtained  as  solutions  of  a  recursive  system  of 
first-order  linear   ordinary  differential  eq\iations.   Initial  conditions  to 
fix  the  solutions  of  this  recursive  system  CEm  be  specified  in  some  cases. 
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I.  Introduction 

In  view  of  the  difficulty  in  olitainlng  exact  solutions  for  problems  in  partial 
differential  eatiations  v;ith  prescribed  initial  and  botindary  conditions,  attention 
has  been  given  in  the  literat-'xre  to  solutions  in  the  form  of  asymptotic  series.  The 
present  paper  presents  a  method  for  obtaining  asymptotic  solutions  for  certain  classes 
of  problems  involving  linear  second-order  hyperbolic  differential  equations.  The 
problems  considered  arise  directly  in  physical  sittiations  or  result  from  reducing 
vectorial  problems  of  electromagnetic  theory  to  scalar  problems. 

Briefly  stated, the  contents  of  this  paper  is  as  follows.  We  consider  the  linear 
hyperbolic  partial  differential  equation 

(1.1)     L(u)  =  ^   a^^i^  +  51  h  Uj^  +  cu  =  0, 
i,j=l      "^   k=l 

il   k 
wherein  u  is  a  function  of  x^ ,  i„,...,x  ;  the  coefficients  a  ",  b  ,  and  c  are  functions 

of  T,  ,  T,.....,x     ,;  and  u.  ,  =   ^- —  ,  Associated  initial  and  boundnry  conditions  are 
1   2      n— 1       ij   ox.o3c. 

specified  later.  We  need  only  note  here  that  one  of  the  boundary  conditions  is  non- 
homogeneous.  We  also  consider  the  case  wherein  a  special  type  of  non-homogeneous  term, 
representing  a  source  or  applied  force,  appears  on  the  right  hand  side.  We  shall  write 
u(x,t)  for  vi(x, ,Zp,...,x  ^,T   ),  letting  x  stand  for  the  first  n-1  variables  and  t  for 

X  . 

n 

Among  solutions  of  equation  (l.l)  we  consider  a  special  solution,  U(x,t),  called 
a  pulse  solution,  which  corresponds  to  a  source  or  impressed  force  whose  time  behavior 
is  the  Heaviside  unit  function.  We  then  introduce  a  Duhamel  principle  which  expresses 
the  solution  u(x,t)  due  to  a  source  or  impressed  force  with  arbitrary  time  behavior 


in  terms  of  en  integral  involving  U(x,t).   If  we  now  restrict  the  arbitrary  time  be- 

-itut 
havior  to  a  harmonic  one,  najnely  e    ,  wherein  «  is  the  circular  frequency  of  the 

impressed  force,  it  is  poaslble  to  show  first  that  with  increased  time  u(z,t)  takes 
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the  form  u(3c)e    ,  The  function  u(j:)  may  be  described  as  the  spatial  part  of  the 

"stesdy  state"  behavior  of  u(x,t).   It  is  possible  next  by  applying  integration  by 

parts  to  the  Diihamel  integral  to  obtain  an  asymptotic  series  for  u(x).  The  series 

runs  as  follows: 

a  a  (iuj)  a 

The  essential  variable  of  the  series  is  l/w  and  the  successive  terms  contain  increasing 
positive  integral  pov/ers  of  l/w.  The  first  term  of  the  series  is  independent  of  «"  and 
hence  is  the  limiting:  value  as  U)  ->oo.  This  term  gives,  therefore,  the  "geometrical 
optics"  limit  of  the  qujBntity  u(x).  The  value  of  this  first  term  at  any  point  x  in 
space,  e.part  from  the  phpse  fpctor,  is  the  Sum  of  the  discontinuities  with  respect  to 
t  in  U(x,t)  at  thst  snme  x.  The  essential  part  of  the  coefficient  of  the  second  term 
of  the  series  is  the  sum  of  the  discontinuities  with  respect  to  t  in  U. (x,t)i  and  so 
on  for  the  higher  coefficients. 

The  very  form  of  the  asymptotic  series  tells  us  th^.t  the  behavior  of  solutions  of 
pprtial  differential  equations  is  Intimately  related  to  the  discontinuities  of  the  pulse 
solution.  However  the  primary  objective  of  this  paper  is  to  show  how  to  obtain  the 
coefficients  of  the  expansion  without  knowing  the  pulse  solution  itself.  For  this 
purpose  we  study  discontinuities  of  solutions  of  eq.ua.tion  (l.l).  To  treat  discontinu- 
ous solutions  we  replace  the  differential  equation  by  a  more  general  integral  equation 
which  admits  discontinuous  solutions,  and  derive  conditions  which  the  discontinuities 
in  U,  U.  ,  U.^,  etc., must  satisfy. 

There  results  first  of  all  the  fact  that  if  such  discontinuities  lie  on  a  surface 
0(x,  ,,..,x  )  =  0  then  fl   must  satisfy  the  characteristic  condition 


(1.3)     Yl       ^   ^Ji   =  0         on  ji  =  0. 
1J=1       ^ 


The  most  useful  conditions  which  follow  from  the  study  of  discontinuities  pre  what 
may  be  called  the  transport  and  higher  transport  equations.   If  one  follows  a  die- 
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continuity  in  IT,  [uj  ,  alon/'  a  path  on  jfl  =  0  which  has  at  each  point  the  direction 

n   . 
callftd  the  transverse  direction,  i.e.,  ^.(s)  =  XI  ^  ^v»  *^®ri  alon^  this  path[ul 

^      k=l    ^ 
satisfies  the  transport  eq[xmtion 

i,,i=l    ^"^   k=l   ^         "^ 
Likewise,  a  discontinuity  in  U.  satisfies  the  hifher  transport  equation 

(1.5)  (  t    a%  .  t.\)[nj.  2  ^=  -  f  a^^  1^-  eV^-  cM. 

i,j-i      ^J     k=i     ^     *^        s^  i7J=i      ^^^^j     k^i      ^^-k        ^ 

15ntirely  analogous  equations  hold  for  the  discontinuities  in  U. . ,  etc.  Solution  of 
these  transport  equations  wo-jld  give  the  values  of  the  discontinuities  required  in 
the  ftsymptotic  series  (1.2).   Initial  conditions  required  to  fix  the  solutions  of 
these  transport  equations  cannot  he  prescribed  generally, hut  some  indications  of  what 
can  he  used  in  special  cases  ere  f^iven. 

A  word  should  he  said  ahout  the  range  of  validity  of  the  asymptotic  series  of  this 
paper.   It  should  he  eacpected  that  a  series  in  positive  integral  powers  of  —  cannot  he 
applicahle  everywhere.   It  cannot  he  expected  to  hold  on  the  houndary  of  shadow  regions 
nor  in  the  neighhorhood  of  foci  and  caustics.   Other  expansions  involving  fractional 
and  negative  powers  of  —  must  ohtain  in  such  regions.   Such  expansions  are  heing  con- 
sidt-redjP.nd  a  forthcoming  paper  hy  Mr.  Irvin  Kay  will  contain  some  results  on  these  ex- 
pansions. Another  point  now  heing  investigated  is  a  mathematical  criterion  for  the 
range  of  validity  of  the  present  expansion.   We  give  here  a  criterion  in  terms 
of  the  properties  of  U(x,t).  That  U(x,t)  possesses  these  properties  may  he  known  on 
jiiysieal  gro^rnds  hut  it  cannot  he  ascertained  at  present  hy  examination  of  the  differ- 
ential equation  and  associated  initial  and  houndary  conditions. 

The  Sijggestion  for  this  work  derives  from  an  analogous  development  carried  out  hy 

the  writer  for  solutions  of  Maxwell's  equations  following  ideas  introduced  hy  R.K. 
Luneherg.  The  present  paper  treats  one  second-order  partial  differential  equation 

1.  See  Kline, M.:  Asymptotic  Solution  of  Maxwell's  Equations .Research  Report  lTo.EM-24 
suhinitted  hy  New  York  University  to  the  Geophysics  Research  Division  of  the  Air  Force 
Camhridge  Research  Lahoratories  or  Comimuiications  on  Pure  and  Applied  Mathematics, 
Vol.  IV,  1951,  pp.  225-262. 
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as  opposed  to  the  system  of  first-order  Mazwell  equations.  While  the  general 
principles  underlyian:  this  development  are  the  sajie  as  those  previously  employed, 
the  details  of  proof  are  quite  different  and  the  resulting  ordinsry  differential 
equations  for  the  coefficients  of  the  ejcpansion  are  also  quite  different  in  form. 
It  is  hoped  that  the  present  material  will  he  easier  to  apply  because  the  ordinary 
differential  equations  to  he  solved  are  simpler.  "For  example ,  esch  ordinary  differ- 
ential equation  here  is  the  analogue  of  t»fo  sets  of  throe  ordinary  differential  equa- 
tions in  the  electromagnetic  case. 

Aside  from  papers  currently  sppearin^  on  the  asymptotic  solution  of  ordinary 
differential  equations  and  on  various  types  of  non-linear  partial  differential 

equations,  several  ha.ve  appeared  recently  which  hear  directly  on  this  paper  or  on 

2 

the  analogous  electromagnetic  casa  .  The  general  haclcground  for  the  exposition  in 

this  paper  will  he  fo^ind  in  Courant-Hilhert :  Methoden  der  Mathem^tischen  Physik 
(Springer,  Berlin,  193?),  V. II,  Chapter  6,  Sections  1  and  2  and  Appendix,  section  k. 


2,  See,  for  example,  Priedlander,  F.G.:  rxeometrical  Optics  and  Maxwell's  Equations, 
Proc.  Cambridge  Phil.  Soc.,  Y,hJ,   Part  2,  Apr.  19^7,  pp.28V286;  Sremmer,  H. : 
The  Jumps  of  Discontinuous  Solutions  of  the  V/ave  Equation,  Communications  on  Pure 
and  Applied  MatheJiatics,  7.1¥,  1951,  pp.  ^19-^26;  Copson,  E.T.:   The  Transport 
of  Discontinuities  in  an  Electromagnetic  Field,  Ihid,  pp. ^27-^33;  Lowell, 
Sherman  C:  The  Propagation  of  Waves  in  Shallow  Water,  Ihid,  V.II,19U9,  pp. 
275-291. 
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II.  Formulation  of  the  Problem 

We  shall  he  concerned  with  ohtPining  asymptotic  expansions  for  functions 
■a(x,  ,  X  ,  •••,j:  )  sfitisfyin^  certain  second-order  partial  differential  equations 
with  associnted  initial  and  hoimdary  conditions.   There  are  two  types  of  problems 
to  which  the  theory  of  asymptotic  soliitions  would  apply. 

Problem  I.   Consider  the  homogeneous  linear  hyperbolic  second-order  partial  differ- 
ential equation 

n  n    , 

(2.1)  l(u)  =  Yl       »  %  1  +  H  D  ^V  +  cu  =  0. 

i,j=l     ^-^   k=l    *^ 

wherein  u  is  a  function  of  3c,  ,x^...,i  ,  and  the  coefficients  are  functions  of 

12  n 

X,  ,   x^,...,x     ,.     The  function  u  is   rea^uired  to  have  all  first-  and  second- 

1  d  n— J. 

order  derivatives  in  some  domain  of  n-dimensional  space.   The  coefficients  a  "^ 
are  likewise  required  to  possess  all  first- and  second-order  derivatives  in 
this  domain. 

Let  S   be  a  hypersurface  in  the  (n-l)-dimensional  space  of  the  variables  z^ , 

x„,«'«,x  ,,  S,   mav  be  infinite  in  extent  or  closed.   In  the  latter  case  we  shall 

2  n-i   1 

be  interested  in  the  solution  u  in  the  space  on  and  exterior  to  SL  .  As  initial 
conditions  we  require  that 

(2.1a)    u(x,t)  =  0  and  u^(x,t)  =  0  for  t  <  0  and  for  all  x  . 

(2.1b)    u(x,0)  =  0  and  u. (x,0)  =  0  f or  x  off  S,  . 

"  1 

djr  boijindary  condition  is 

(2.1c)  u(x,t)   =  h(x)f(t)   for  x  on  S^ . 

It  is  assumed  that  f (t)  =  0  for  t  <  0  and  that  f(t)  and  f '(t)  are  continuous. 

For  this  problem  I  we  define  as  the  pulse  solution  U(x,t)  that  solution  which 
satisfies  the  conditions  (2.1),  (2.1a),  (2.1b), and  (2.1c)  except  that  f(t)  in  (2.1c) 
is  the  Heaviside  unit  function  Vj  (t)^w}iich  is  0  for  t  <  0  and  1  for  t  >  0.  We  re- 
quire further  of  U(x,t)  that  every  bounded  domain  of  (x,t)  space  be  decomposable  into 
a  finite  number  of  closed  subdomains  in  each  of  which  U,  U. ,  U, . ,  and  the  higher  time- 
derivatives  exist  and  possess  ell  first-  and  second-order  continuous  derivatives.   More- 
over all  discontinuities  in  IT,  U.  ,  IJ..  .etc., with  respect  to  t  must  be  finite. 
Problem  II.   Consider  the  non-homog'eneous  linear,  hyperbolic  second-order  partial 
differential  equation 

(2.2)  L(u)  =  h(x)f(t), 

where  h(x)  is  not  zero  in  some  closed  (n-l)-dimensional  region  V  of  x^ ,  • • • ,  x  - 

1        n— 1 

space,  f(t)  =  0  for  t  <  0  and  f(t)  and  f '(t)  are  continuous  for  t  >  0.  The  initial 
conditions  are: 


-  6  - 

(2.2a)   u(T,t)  =  0  and  \\^(x,t)  =  0  for  t  <  0  and  for  all  x  . 

(2.2b)   u(x,0)  =  0  end  u^(x,0)  =  0  for  x   outside  of  V. 

We  may,  if  we  wish,  impose  the  following  "boundary  condition:   Let  S   be  any 
surfece  in  (n-l)-dimensional  space: 

(2.2c)      u(x,t)  =  0  for  x   on  S  and  for  all  t. 

Tot   problem  II  we  also  define  as  the  pulse  solution  the  function  lT(x,t)  which 
satisfies  (2.2)  except  that  f(t)  is  the  specie!  function  "n  (t),  and  U(x,t)  satisfies 
the  initiPl  and  boundrry  conditions  (2.2a),  (2.2b),  and  (2.2c).   The  additional  con- 
ditions imposed  on  U(x,t)  in  problem  I  ps  to  the  existence  of  derivatives  in  closed 
subdomains  and  as  to  the  finiteness  of  the  discontinuities  in  t  are  imposed  here  too. 
Problem  III.   This  problem  is  essentially  a  combination  of  the  two  preceding  ones. 
We  consider  the  non-homogeneous  linear  hyperbolic  second-order  partial  differen- 
tial equation 

(2.3)  ^M   =  h^(3c)f;L(t)  . 

where  h- (x)  ^  0  in  some  closed  (n-l)-dimensional  volume  V  of  ^]^""»^n-l  ^P^°® 

bounded  by  the  surface  S,  f^(t)  =  0  for  t  <  0,  and  f (t)  and  f '(t)  continuous  for 
t  >  0.   The  initial  conditions  are 

(2."^)   u(x,t)  =  0  and  u^(x,o)  =  0  f or  t  <  0  and  for  all  x; 

(2.''b)   u(x,0)  =  0  end  u.  (x,0)  =  0  for  x   outside  of  V. 

In  addition,  let  S   be  any  closed  surface  on  which  we  impose  the  boundary  con- 
dition 


(2.'^c)   u(x,t)  =  h„(x)f„(t)  for  x 


on 


1' 


where  f-(t)  end  f  '(t)  are  continuous  and  f_(t)  =  0  for  t  <  O  . 

Problem  III  can  be  reduced  to  i^roblems  I  and  II,  We  have  but  to  let 


u  =  u,+  u„  where 


L(u^)  =  h^(x)f^(t), 


5u 


cu-i 
Uj^(x,t)  =  -^  (x,t)  =  0  for  t  <  0  and  for  all  x, 


3ui 


u,(x,0)  =  — rr  (x,0)  =  0  for  x  outside  of  V, 
u-(x,t)  =  0  for  X  on  S  and  for  all  t; 


?'nd  •'^^^2^  "  °' 


3u 

n„(x,t)  =  -rr  (x,t)  =  0  for  t  <  0  and  for  all  x, 
c.  ot 

u,(x,0)  =  -rf  (x,0)  =  0  for  X  outside  of  ?, 

U2(T,t)  =  h  (x)  f  (t)  for  X  on  S^. 
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III.  The  Duhamel  Principle 

The  pulsfi  solution  we  have  Introduced  for  the  several  prohlems  will  generally 
he  discontinuous  nt  some  values  of  x  and  t  because  it  corresponds  to  a  sudden  dis- 
turhance  introduced  at  t  =  0  either  on  some  boundary  S   (problem  I)  or  in  some 

region  T  (problem  II).  This  disturbance  will  propagate  as  a  discontinuity  which 
will  reach  point  x  off  S  ,  or  outside  of  V,  at  some  later  time  t.  We  have  re- 
quired that  in  any  bounded  region  of  (x,t)  space  there  be  closed  subdomains  within 
which  U  and  its  derivatives  to  second  order  be  continuous.  However  across  the 
boundaries  of  such  closed  subdomains,  U,  U. ,  U. . ,  etc.,  may  be  discontinuous.  We 

shall  suppose  that  the  discontinuities  lie  on  some  hypersurface  0(x- ,X-, . . . ,x  ^,t)  =0, 
that  0=0  may  be  solved  explicitly  for  t  In  the  neighborhood  of  a  given  value  t^  of 
t,  end  that  the  explicit  solution  is  given  by  t^  =  v/(x^...,x  ^),  where  \ii   possesses 
first-  and  second-order  derivPtives  with  respect  to  its  variables. 

Since  U  is  discontinuous  across  the  boundaries  of  the  closed  siibdomains  we  shall 
require  that  certain  discontinu: ty  conditions,  (3.11)  and  (3.12),  hold  on  jt>  =  0. 
Under  these  conditions  and  the  general  assximptions  made  about  U(x,t)j  the  following 
assertion  holds: 


(3.1)      u(x,t)  =  l^y  u(x,t-r)f('r)  dT 


This  integral  which  relptes  the  solutions  u(x,t)  and  U(x,t)  is  known  as  Duhamel's 
integral". 

We  shall  verify  that  (3.1)  is  correct  for  problem  I  and  then  note  any  special 
considerations  that  me-y  be  needed  to  take  care  of  problem  II. 

In  verifying  (3.1)  we  shall  take  into  account  a  discontinuity  at  t  =  t^  of  U(x,t). 
let  us  suppose,  then,  that  U(r,t)  is  discontinuous  at  t  =  t-  and  that  0  <  t^  <  t. 
Since  we  have  required  that  TT(x,t)  be  sectionally  smooth  and  that  it  satisfy  (2.1)  in 
each  closed  domain  in  which  it  is  sectiona-lly  smooth, we  shall  substitute  (3.1)  into 
(2.1).  "First,  we  note  that  from  (3-1)  we  may  write 


(3.2)        f     u(z,s)ds  =  J^   U(x,t--y)f(T)d'T'  . 


0  "0 

In  (2.1)  we  write  s  for  t  and  integrate  with  respect  to  s  from  0  to  t.   Since  the 
coefficients  of  (2.1)  are  assumed  to  be  independent  of  t,  we  may  write 


3.  Of.  CoTirant,  R.  and  Hilbert,  D:  Methoden  der  Mathematiachen  Physik,  V.II,  p. 18^, 
Springer,  Berlin,  1937. 
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;3.3)  ^  a^-'  ^-r—    /   u(x,g)ds  +  ^  IJ  rf-  I  u(x,s)ds  +  c  /  u(x.s)ds  =  0. 


We  may  now  use  (3.2)  to  write 
n     ..   .2    ^t 


(3.^)  H   a-«5  —^  f   TT(x.t-1')f('V)dr+  ^b'^  -J-  [  uU.t-T)f(r)dT 

1  u(x,t-r)f('Y)d'Y  =  0. 


+  c 

0 


Since  U(x,t)  is  discontinuous  at  t  =  t^  ,  lT(x,t-t')  is  discontinuoiis  at  t-T'=  t^  or 

1*  =  t-t,  •  V/e  therefore  shall  treak  up  the  Y  -  interval  of  integration!  into  0  to 

t-t^-  and  t-t^+  to  t,  where  by  t-t^-  we  mean  t-t^-e  ana  by  t^-t^+  we  mean  t-t-,+€. 

We  shfill  use  the  fact  th^^t  t,  =  wCx^  ,x„,  ♦  •  •  ,x  ,), 

1      12      n-1 

To  verify  that  (3.^')  is  correct  we  shall  consider  the  effe^^t  of  the  differen- 
tiations in  each  term  of  (3.^).   Consider  first  a  term  such  as 

ii      a^       .t 

""     sx  a-    J  u(^.^^-'>Of(r)dr,  for  ij  =  1.2,- ••,  n-l. 

i  "J       0 


]    U(x,t-' 


^'^'  Sjt  J  n(x.t-r)f(r)dr 

1     J       o 


-t-t,- 

=  a 


^^  r-^-  r      ^  U(x,t-T)f  (T)dT+  a^''  — ^  f      UCx.t-Tlf  (T^dT 

ax.?vx^^  ^^i^^j^t-t^ 

r ;-« - 

-  \  J       Uj(x,t-  T)f  (T)dr-U(x,t^+)f  (t-t^)  vj 

+y      Uj(x,t-'r)f(r)dT+u(T,t^-)f(t-t^)  V 


/  . 


We  now  carry  out  the  differentiation  with  respect  to  x . .   In  differentiating  the 

second  and  fourth  terms  we  note  that  we  have  a  product  of  three  f-unctions  in  which 
product  the  first  term  is  a  compound  function  of  x. .  Hence 
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0  4-  4-      +■  — 

^■^  ^^TtJ  n(x,t-r)f(Y)dr=  a^-'J     \j(x,t-T)f(f  )dr -Uj(x,t^+)f(t-t^)x„. 


t-tl+ 


-U(x.t,  +  )f(t-tjM/,  .  +   f  TJ      (3r,t-T)f(Y)d'Y+  U.  (x,  t.-)f  (t-t  J^u  . 

1  i     1,1      ^;  t-ti+     ♦^  J         -L  -L     1 

+n(x,t^-)f(t-t^>i/. . 


V?e  may  shorten  this   result   somewhat  hy  using   the  notation   [TT(x,t^  ]]    to  mean  U(x,t^  +  ) 
-  U(x,t^-).      Then,   for   i,j   =   l,2,...,n-l, 

2    /-t  •  •[  r'*^ 

(3.5)  ^  \^\^      J    U(x.t-Y)f(r)dY=  a^^l  j  U.j(x,t-'r)dT+[y(x,t^}]f' (t-t^)v.,,,^ 


-  Iu(x,t^j]f(t-t^)M/.j  -  [U^(x,t^[lf(t-t^>l/j  -  Il\(x,t^2  f(t-t^Vi 

-  [n^(x,t^[l  f(t-t^)v^wj 

¥e  consider  next  terms  of  the  first  sum  in  (3.^)  for  which  J  =  n  and  i  Is 
1,2, . . . ,n-l. 

"t  .2   r^-H-  .,  ,2  „t 

0 


''^/o^^(-'^-'^^^('^)'^'^=  ^''ial/o  ^J(-.t-r)f(r)dr.  a^^^/  u(x.t-Y)f  ('r)dr  , 


i   "  t-t^+ 


=  a^*-^(      ujx,t-r)f(r)dr4  u(x.t,+)f(t-t-)+  u.(x,t-r)f(Y)  dr 

^»i  r  0     *  1  ^    J  t-t^+  ^ 

+TT(x,0)f(t)  -  Tl(x,tj^-)f (t-t^) 

Before  cprryin^  out  the  next  differentirition  we  shpll  drop  the  term  containiAP:  U(x,o), 
since  "by  the  initial  conditions  this  is  0  except  when  x  is  on  hoijndary  surface  S  .  The 
latter  case  will  he  considered  separately. 


Now 


2  A 


9x.9t  % 


dr  =  a 


it 
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,t-t,- 


V.  , 


t-t  + 


-U^(x,t^-)f(t-t^)M/^+U(x,t^-)f' (t-t^V. 


(3.6)  a^*^^  f  ij(x,t-y)f(r)dr  =  a^*  ■  f  u.^(x,t-'r)f('>')d'v-[u(x.t^2f '(t-t^)v^ 

Those  second-order  terms  of  (3.^)  for  which  i  =  n  and  J  =  l,2,...,n-l  g;ive  the 
same  result  as  ( 3. 6) .since  a  "^  =  a''  ,  and  the  order  of  the  differentiation  is  immaterial  , 

We  now  consider  the  second-order  term  in  (3.^)  for  which  i  =  J  =  n.  The  first 
of  the  differentiations  to  he  performed  in  this  connection  was  done  in  deriving  (3.6), 
Hence 


.2  „t 


t-U- 


a 


**-^  J  TT(x.t-r)f  (r)dr  =  a"|^  I  j     ^  ujx.t-r)f  (V)  dt  +u(x,t^+)f  (t-t^) 

5t    0  o 


f        u,  (x,t-'r)f('v)dr-u(x,tT-)f(t-t  ) 
•^  t-t,+  *  ^ 


tt, 

=  a 


(.t-t,- 


+  J         U^^(x,t-r)f('r)dr+U^(x,0)f(t)-U^(x,t^-)f(t-t^)-IJ(x,t^-)f'(t-t^) 

We  use  the  initial  condition  that  U. (x,0)  =  0  for  x  off  the  boundary  surface  S^,and 

oh  ta  i  n 

U.  In  performing  the  differentiation  with  respect  to  x.  first  and  then  with  respect 
to  t  one  does  end  up  with  the  term  U.(x,0)  f(t)in  addition  to  the  other  terms  in 
(3.6).  However  since  'J(x,t)  =  0  for  t  =  0  except  on  S,,  we  conclude  that 
Uj(x,0)  =  0  . 
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tt  a' 

a    ; 

at' 
(3.-?) 


J  u(x.t-r)f(V)d'r'=  a**  I  J  U^^(x,t-y)f('r)dY+  [U(x.t^)]f' (t-t^) 

*  [nt(x,ti)]f(t-t^)  r . 


To  perform  the  differentia tiocs  called  for  in  the  first-order  terms  in  (3.^)  we 
may  use  the  results  already  ohtained.  For  k  =  l,2,...,n-l, 

(3.8)    h^^  jTiu,t-r)f(Y)d'T'  =  •b'^l  J  u^Cx.t-iOf  (Y)d'y  -  [u(^.til]f  (t-t^vj  i  . 


For  k  =  n  v/e  have 


(3.9)  b*^  f  n(x,t-V)f(Y)dy=  h*   J  U^U,t-'Y)f(V)dV+  [U(x,t^i]f(t-t^)  I  . 

I  °  J 

Tor  the  remsininf  term  in  (3.^)  we  merely  restate  that 

(3.10)     c  J  u(x,t-r)f(r)dr  =  c  I"  uCx.t-T^f  (r)dr  . 

0  0 

To  verify  the  correctness  of  the  Duhamel  principle,  that  is,  that  (3.^^)  is 
correct,  we  must  now  pdd  up  the  terms  in  equations  (3.5)  to  (3.10).  We  find  on 
the  right  side  th8.t  the  terms  involving  integrals  amount  to 

J     L(u(x,t-'0)  f(r)dV, 

Since  U(j',t)  is  required  to  satisfy  L(u)  =  0  for  t-Y<  t^  and  t-Y>  t^  ,  these  terms 
add  up  to  zero. 

V/e  collect  next  those  terms  on  the  right  side  of  equations  (3.5)  to  (3.10)  in- 


volving f '(t-t.).  These  are 


f  (t-tj^)^ 


n-1   .  , 


n-1 


n-1 


XI  a^''[ll(x.t-llHr.u/.  -  IIa"[u(x,t  3^4,  -  H  a  "^  &(x.tT  2^^ 
i,,i=l        ■'^•'1=1         ^   ^   j=i        ^   J 

+  a"lu(^.tlll[   . 


We  recall  that  t  =  \vix^,..,x     ,)  is  an  explicit  form  of  0(Xt,x„...,z  )  =  0.   Hence 
J  ^  ^   \^  n-1  '^  1  2     n 

\y.  =  -  ■7—  .  We  see  therefore  th-'it  these  terms  amount  to 
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f'Ct-tJ  n     ,. 

0^  ij=l      -^ 


We  require  that 
n 


(3.11)    Yl    ^^   0.  0.  =  0 
i.j=l       -^ 


i^A  d(  -f^on^  =  0; 


thnt  is,  if  0  =  0  is  a  surface  on  which  U(x,t)  is  discontinuous,  then  0  is  required 

to  satisfy  this  first-order  differential  condition  on  0  =  0.  This  condition  is  known 

as  the  characteristic  condition  . 

\Ie   collect  next  the  terms  involving  f(t-t^).   Omitting  this  factor  these  terms 

are:  ^  .  ^ 

n-1      .  .  n-1      ,  .  n-1      .  . 

-  H     a^-'ljKx.t    l]v       -  ^    a^J[u.(x,t    [Iv/.   -  E    a^-'lu   (z,t J]iv. 
i,J=l  ''     i.j=l  -"      i.,i=l  "^ 

-E  a^^[u.(x,t.aM'.v, +  ra"[:u.(x.t.i]  +  z:a*''[u  (x.t.a 

i,j=l  ^^       i=l  j=l  J         -^ 

+  a"[:U^(y.t^a    -  ZI  ^''&(x.t^))^^  +  b^[u(x,t^[l       . 


We  shall  show  later,  in  connection  with  the  first  discontinuity  condition  of  section 
V,  that  these  terms  can  he  rev/ritten  as  on  the  left  side  of  (3.12).  We  now  require 
that 

(3.12)  E  C^]  a'''0ii  +  t  Cu:ih\  +  2  JZ    a'^Qjri  0  =  0  on  0  =  0. 
i,j=l        ^"^   k=l        "^    i,j=l         -^ 

on 

Conditions  (3.11)  and  (3.12)  pre  imposed  on  U(x,t)  and^he  discontinuity  surface  0  =0 
in  order  thnt  the  Buhamel  principle  hold  . 

V/e  must  now  verify  that  the  Duhamel  principle  holds  for  points  x  on  the  boundary 
surface  S^.  On  this  surface  u(x,t)  =  h(x)f(t),  whereas  Tj(x,t)  =  h(x)>i(t).  Then 

h(x)f(t)  =  1^  J^h{x)>^(t-r)fi'v)  dry. 

Since  ^[(t-'rt  =  1  for  0<'V<  t,  the  verifier tion  is  immediate   not  only  for  t>0  hut  as 
t  -^  0+.   Likewise  the  initial  conditions  are  satisfied. 

5.  See  Courant,R.and  Hilhert,D.,  loc.cit. ,V. II,p, 351.  The  significance  of  (3.11)  and. 
(3.12)  will  he  cler^rer  after  the  discussion  of  section  7  of  this  paper. See  the  dis- 
cussion of  equation  (5.12)  and  see  (5.15). 

6.  See  the  preceding  footnote. 
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We  consider  next  the  modifications  needed  in  the  pljove   verification  of  the 
Duhpmel  principle  to  apply  it  to  Problem  II.  Equations  (3.5)  to  (3.10)  are  still 
cori'ect  statements.   Moreover  we  shall  again  impose  conditions  (3.11)  and  (3.12), 
so  that  those  terns  which  were  previously  discarded  on  the  right  sides  are  a^ain 
discarded.  However,  the  stim  of  those  terms  on  the  left  sides  of  equations  (3.5) 
to  (3.10)  which  involve  an  integral  now  yields,  in  view  of  the  non-homogeneous  term 
In  {?,.?)   and  step  (3.3), 


/: 


it 

h(x)f(3)d9. 


o 

The  terms  on  the  right  side  of  steps  (3.5)  to  (3.10)  which  contain  the  integral 
of  TJ  and  its  derivatives  amount  to 

(3.13)        J    L^J(x,t-7'))f(T)d'r  . 

Since  U  is  required  to  satisfy  the  equation  L(U)  =  h(x)>»(t),  the  quantity  (3.13) 
equals        . 

Since  •h(t-T)  is  1  for  0  <  'X<  t,  we  see  that  the  Duhsmel  integral  for  u(x,t)  sat- 
isfies the  differential  equation. 

As  to  the  initial  and  houndary  conditions,  satisfaction  of  the  differential 
ea.ua tion  takes  care  of  points  x  not  only  outside  of  7  hut  also  in  7,  for  h(x)  /  0 
only  in  V.  Given,  also,  a  houndnry  condition  that  u(x,t)  =  0  on  some  surface  S^^, 
we  require  that  'T(x,t)  satisfy  the  same  condition.   It  follows  at  once  from  (3-1) 
that  for  X  on  S^ ,  the  left  and  right  sides  are  both  zero  for  all  t  >  0  and  es   t 
approaches  0+. 


_  IL  _ 

17.  Derivation  of  the  Asymptotic  Expansion 

We  now  restrict  the  time  behnvior  of  the  f (t)  which  appears  in  the  houndrry 

— iwt 
condition  of  Prohlem  I  and  the  differential  equation  of  Prohlem  II  to  e     and 

we  -use  the  IKitejiiel  integral  to  derive  the  asymptotic  expansion  for  u(x,t).  Strictly 

we  shsll  show  that  u(y,t)  contains  p  "transient"  term  which  approaches  zero  rs   t 

"becomes  infinite  and  that  the  remaining  "steady  state"  pert  of  u(x,t)  is  of  the  form 

u(x^,x_, • •• ,x  _, )e    .   Ve  then  derive  an  asymptotic  series  for  u. 

We  begin  with  the  Dohamel  integral^ where in  we  let  f (t)  be  e~   .  Then 
(it.l)         u(x,t)  =  1^  /'u(x,t-'r)e"^"'V  . 
Let  t-'V  =  s  and  t  =  t.  With  this  change  of  variable 

a  C^„f     X  -iu)(t-s)^ 

Ti  =  -^J    U(x,8)e       'ds  , 

0 

~  at  ^®    J   ^(^-•^^®   ^^^  • 

0 

=  U(x,t)-lo)e"^"'*  y  U(x,s)e^'"''ds. 

Adding  and  subtracting  equivalent  terms  gives 

-t 


u  =  U(x,t)-U(x,  oo)-iu'e"^'"*  /"   (u(x,s)-U(x,  co) J  e^"'^ds+U(x,  oo) 

U(x,oo)-iw   I    (u(x,s)-U(x,  oo)j 


-iUJt 

e 


=  U(x,t)-lT(x,oo)+e~^"'*  1 


iOJS^ 

e       ds 


We  now  let  t  become  infinite.  We  shell  neglect  whpt  may  be  considered  on  physical 
grounds  as  a  transient,  namely  the  terms  Tj(y,t)-Il(x,  oo),  because  at  any  point  x,  the 
field  due  to  the  pulse  should  esta>>lish  itself  in  a  "short"  time  and  readily  approach 
a  static  condition.  ?or  the  same  rihysical  reason  too  we  expect  that  the  integrand 
should  approach  zero  rapidly  for  large  b  and  hence  that  the  Integral  from  0  to  oo 
should  not  only  exist  but  differ  little  from  the  integral  from  0  to  t.  Hence  with 
increasing  timejU(x,t)  taices  the  form  u(x)e    ,  and 

u  =  U(x,  go)-!"'  f     (u(x,r)-U(x,  od))  e^^  dY  . 

We  now  integrate  by  parts  to  obtain  the  asymptotic  expansion  for  u.   One  such 
Integration  yields 
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CX)    /lOO 


u  =  LT(x,oo)  -  (u(x,'V)-U(x.oon  e   I  +yQ  U^(x,Y)e   dt  . 

In  evaluating-  the  second  term  we  must  take  into  eccoimt  the  discontinuities  in 
U(x,t)  at  values  of  t.   We  have  assumed  that  these  discontinuities  are  finite. 
There  may  be  an  infinite  number  of  them,   let  Y  ,  Y  ,  •••,  'V  ,•••  be  these  dis- 
continuities  and  let  lu(x,'V  []  =  LT(x,T +)-U(3!:,Y -).  Then 


a*-"        a       a 


a  0 

Another  inteprstlon  'hy  parts  yields 

a  a.  *'o 

The  second  snmitatlon  Includes  U^  (x,0)  and  U^  (x,oo).  We  may  expect  on  physical 
grounds  that  U  (x,oo)  is  zero,  because  U(x,t)  approaches  a  static  field  as  t  in- 
creases. The  quantity  U  (x,0)  should  also  be  zero  at  points  x  away  from  the  source, 
since  tJ(x,t)  is  0  there  for  t  less  than  some  positive  constant.  Hence  we  (?i5(card 

these  values  of  t  in  the  second  summation. 

rr 

It  is  apjjarent  that  the  process  of  integration  by  points  can  be  continued  and 
the.t  we  obtain  a  series  in  positive  integral  powers  of  —  for  u.  We  therefore  write 

iujV  ittJY  iiuY 

a  a  ec 

For  points  x  not  on  S.  or  at  the  source,  U(x,0)  equals  0  and  may  lie  neglected. 

This  infinite  series  is  certainly  asj-mptotic  to  u  since  the  integral  neglected 
after  any  stape  is  seen  by  an  integration  by  parts  to  be  of  degree  one  higher  in  — 
than  the  preceding  ones.   The  argument  assumes  tha.t  the  integral  converges.  This 
may  be  expected  on  physical  grounds  ,  because  the  higher  derivatives  of  U(x,t)  should 
certainly  approach  zero  rapidly  with  increasing  t.   If  we  stop  the  integration  by 
parts  at  any  one  term  we  obtain  an  exact  esq^ression  for  u(x)  wherein  the  remainder, 
that  is,  the  integral, is  of  order  one  higher  in  —  than  the  last  term  of  the  finite 
series. 

At  an.y  x,  the  f  in  (^.2)  are  the  values  of  t  at  which  U(x,t)  or  any  one  of  its 
time-derivatives  is  discontinuous.  We  have  already  noted  in  connection  with  t^e  veri- 
fication of  the  Duhamel  integral  that  if  (x,t)  is  a  point  in  the  n-dimensional  space 


7.  The  argument  here  shoiild  later  be  replaced  by  a  ma.thematical  one  v.'hich  wotild  enable 
us  to  decide  from  a  knowledge  of  the  differential  equation  itself  and  the  initial 
and  boimdary  conditions  whether  the  discontinuities  in  the  time  derivatives  of  lT(x,t) 
remain  finite  and  that  these  time  derivatives  approach  zero  rapidly. 
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at  which  TT(7:,t)  is  discontinuous  .then  (x,t)  lies  on  a  hypersurface  0  =  0  which  we  h^ve 
required  to  sati'^fy  the  characteristic  condition  (3.11).   It  will  be  shown  later  thR.t 
the  points  (x,t)  at  which  U.  or  any  hi,eher  time-derivative  is  discontinuous  likewise 
lip  on  (i  =  0.  Fence  ^iven  an  x,  the  corresponding  'V,  inay  l^e  determined  hy  solving 
f)   =  0  for  t.  There  may  of  course  he  many  T^  for  a  given  x  and  we  shall  therefore 
denote  the  explicit  solutions  hy  T^  =  \j/  (x,,x„,'»*,x  ,).  Hence  we  may  replace 
the  ^^   in  (^.2)  by  the  corresponding  x^   , 

It  is  obvious  from  the  form  of  (^.2)  that  the  asymi-totic  series  should  be  useful 
for  large  uj  and  that  for  w  =  oo,  the  first  term  yields  the  geometrical  optics  field 
in  the  case  of  electromagnetic  waves  and  the  corresponding  limits  in  the  cases  of  other 
physical  phenomena.  Actually  the  series  (^.2)  gives  more  than  the  usual  geometrical 
optics  methods  do,  for  the  phase  factors  e   °'   are  also  included. 

To  use  the  series  (^.2)  one  must  know  the  coefficients.  We  have  already  indicated 
that  we  know  the  V  .  However  the  essential  part  of  each  coefficient  is  a   discontinuity 
or  sum  of  discontliuiiti es  of  U(x,t)  or  one  of  its  time-derivr. lives .  These  are  generally 
not  known  at  various  points  (x,t):  however  it  niay  be  possible  to  find  them  without  re- 
quiring a  knowledge  of  the  entire  solution  U(x,t).  We  shall  therefore  undertake  a 
general  study  of  discontinuities  of  solutions  of  equations  (2.1)  and  (2.2),  from  which 
we  shall  obtain  ordinary  differential  equations  for  each  of  these  discontinuities.  These 
differential  equations  may  then  be  used  to  determine  the  discontinuities  directly  and 
thereby  the  asymptotic  expansion  itself. 


8.  See  section  VI. 
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T.  rerivation  of  the  First  Discontln-uit:-  Condition  and  First  Transport  Eqiiation 

We  proceed  to  study  discontinuities  of  solutions  of  the  differential  equations 
(2.1)  end  (2.2).   Of  course  at  a  point  or  elon^  e>   hypersurface  on  which  tt  or  its 
first  or  second  derivp-tives  are  discontinuous  the  differential  equations  have  no 
mea.ninf'.   Moreover^  if  some  discontinuity  surface  separates  a  ref'ion  G  in  which  (2,1) 
holds  from  a  region  G^  in  which  (2.1)  holds. there  is  no  essential  connection  between 
the  two  solutions  in  the  two  regions.   Initial  end  boundary  conditions  may  determine 
the  solution  of  (2.1)  uniquely  in  ref'ion  G^  but  the  solution  in  G  may  still  be  in- 
dependent of  the  one  in  G  ,  Evidently  to  discuss  discontinuous  solutions  of  these 
differential  equations,  some  relationship  must  be  imposed  to  connect  the  solutions 
on  either  side  of  a  di?!Continuity  surface.  We  shall  discuss  equation  (2.1)  and  rema.ric 
later  on  the  minor  additional  featiu-es  required  to  treat  (2.2). 

One  can  roquire  that  a  discontinuous  solution  of  (2.1)  be  a  limit  of  continuous 

solutions  each  of  which  satisfies  (2.1)  in  both  regions  and  on  the  discontinuity  sur- 

9 
face  .  The  discontinuous  solution  may  then  be  expressed  in  terms  of  the  Hesviside 

unit  fxmctlon  and  substituted  in  the  aifferential  equation.  Under  this  procedure 
derivatives  of  the  unit  function  are  the  6-f unction  and  its  derivatives.  By  suitable 
operation  with  this  function  it  is  possible  to  obtsin  the  discontinuity  condition  ob- 
tained below  and  even  to  secure  the  higher  discontinuity  conditions  obtained  in  the 
next  section  .  The  use  of  the  6-function  and  its  derivatives  can  be  made  n^athemati- 
cally  rigorous  by  using  the  theory  of  distributions  due  principally  to  Laurent  Schwartz  . 

To  achieve  rigor  and  to  avoid  dependence  upon  the  somewhat  elaborate  theory  of 
distributiciB  we  shall  employ  another  procedxire.  We  shall  replace  the  differential 
equation  (2.1)  by  an  integral  eq-uation  which  Is  equivalent  to  (2.1)  in  a  region  G 
in  which  the  first  and  second  derivatives  of  u  exist.  The  integral  will  be  meaning- 
ful, however,  even  for  discontinuous  u, and  we  shall  require  that  the  discontinuous 
solutions  u  considered  in  this  paper  satisfy  the  integral  equation.  From  this  equation 
we  shpll  derive  discontinuity  conditions  which  must  hold  on  any  discontinuity  surface. 
These  discontinuity  conditions  will  then  be  converted  into  ordinary  differential  equa- 
tions or  transport  equations  for  the  dlRContintuities  in  u,  u.  ,  'li^+t  etc.   The  solutions 
of  these  differential  equations  corresponding  to  suitable  initial  conditions  wre  the 
coefficients  of  the  ps^rmptotic  series  (^.2).   In  this  section  we  shall  obtain  the  first 
discontinuity  condition  and  first  transport  ea_uation. 


9.  See  Courant-Hilbert,  loc.  cit.,  p.  36O. 

10.  For  the  case  of  the  vmve  equation,  the  firrt  discontinuity  condition  is  obtained 
by  the  procedure  just  described  in  Bremmer,  H. ,  loc.  clt.   See  the  last  part  of 
this  section. 

11.  Schwartz,  L. ,  Theorie  des  distributions,  2   vols.,  Hermann  and  Co.,  Paris,  I95I-52. 
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Let  u  and  Ite  first  and  second  derive tives 
exist  and  "be  "bounded  in  a  regipr.  G,   Let  IL  "be 
any  prbitrary  function  of  x-,  ,Xp,f  •  ,x  ,x  =t,  such 
that  -TL  and  all  first-order  derivatives  vanish 
on  the  ■boundary  and  such  that  -TL  possesses  de- 
rivatives of  any  finite  order  with  respect  to 

the  T^    in  G.  Let  dv=dx^  dx. .  • • •  dx  .  Then 
1  I       c  n 


(5 


.1)  /  iLL(u)dv  =  /  il  H  a^V  ,dv  +  /  A  E  ^\^^  +  ./-H. 
G  G   i,j      *'  k  "g 


cudv. 


Let  (dx.)  indicate  that  partial  integration  with  respect  to  dx,  has  "been  performed, 

J  J 

Then  by  integration  "by  parts 


Gi,j  ri,j  Gi,j 


dv  . 


!^ince  -XU  vanishes   on     T     ,    the  surface   integral  vanishes.      Another   integration  hy 
parts  plus  the  fact  that  -TL  and  i\,  vanish  on     P    yields 


(5.2) 

Likewise 
(5.?) 


dv. 


/fuEa'\,dv=      /E(-^ai^)     udv. 
G        i.j  ^-  '^G  i.j  ^"^ 

/jXE  ^\dv  =  -  /    E(Ah^)  udv. 
•^G       k  ^  -^  G  k  ^ 

We  may  therefore  state  from    (5.2)  and   (5.3)   that 

(5.^)  yjT.L(u)dv  =   f  H  u(/L.a^'5)^jdv  -   /  E  u(Ab^)j^cLv  +  y  Acu 

G  Gi,j  Gk  G 

12 
Now  the  adjoint  differential  equation  of  L(u)  is,  by   definition, 

M(n-)  =  E  (^a"'').,  -  E(n-b^)v  +  c  a  . 

i.j      ^   k      *" 
Hence  (5.^)  tells  us  that 

(^•^)  /  XLL(u)dv  =  f  uM(rt)dv. 

^G  "^G 

We  conclude  also  from  the  above  steps  that  in  a  region  Gvikerein  u  has  all  first 

and  second  derivntives  which  are  bounded.the  differential  equation  L(u)=0  may  be 

1'^ 
converted  into  the  integral  equation  ' 


(5.6) 


y  u^^(A) 


dv  =  0, 


12.  Courant-Hilbert,   loc.cit.,  p.^3^  . 

13.  ef.    Courant-Hilbert,   loc.cit.,   pp.   i»69-'^0. 
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wherein  it  is  -Luiderstood  thet  u  must  satisfy  the  integral  equetion  for  any   jf).  which 
possesses  first  and  second  derivatlyes  in  &  and  is  such  that  XI   and  its  first  deriv- 
atives vanish  on  the  "boundary  of  G.  Moreover,  since  the  steps  leading  from  (5.1)  to 

(5.5)  pre   reversihle  hy  virtue  of  the  arbitrariness  of -O.,  the  differential  equation 
L(u)  =  0  is  entirely  equivalent  to  the  integral  eq\ia,tion  (5.6). 

If  G  is  a  region  in  which  u  or  any  of  its  first  or  second  derivatives  are  dis- 
continuous, then  the  equation  L(u)  =  0  hps  no  meaning  throughout  G.   However,  equation 

(5.6)  dc^Sj pnd  we  therefore  pgree  to  replace  the  differential  equation  hy  the  more 
general  integral  equation  (5.6).  We  note  that  the  iise  of  (5.6)  as  our  starting  point 
implies  that  L(u)  =  0  holds  in  any  suhregion  of  G  in  which  u  and  its  first  and  second 
derivatives  exist,  for  we  may  take  a  hypersphere  about  any  point  in  this  subregion 
choose  an  iX  which  with  its  first  derivatives  vanishes  on  this  hypersphere  and 
apply  the  steps  leading  from  (5.6)  to  (5.1).  We  require  not  only  that  a  discon- 
tinuous u  satisfy  the  integral  equation  (5.6)  tut  that  u  and  its  first  and  second 
derivatives  be  contimious  from  within  and  bounded  over  those  closed  subregions  in 
which  u  and  its  derivatives  exist.  Otherwise  stated,  u  and  its  first  and  second 
derivatives  must  be  continuous  except  on  a  surface  0(x^,...,x  )  =  0  on  which  all 
discontinuities  of  u  and  its  successive  time-derivatives  may  have  only  jujnp  dis- 
continuities.  It  is  assumed  also  that  the  first  and  second  derivatives  of  0  exist. 

The  conditions  on  u  and  on  Si"   specified  above  suffice  to  derive  the  discontinuity 
conditions  and  first  transport  equation  in  this  section.   However,  to  derive  later 
the  discontinuity  conditions  and  transport  equations  for  [u+J  ,  01^.^]  ♦  etc.,  we  require 
the  further  conditions  that  within  any  bounded  region  G,  there  exist  a  finite  number 
of  closed  subdomains  within  which  u .  ,  u^o..  etc.  and  their  first  and  second  derivatives 
with  respect  to  all  variables  be  continuous  and  bounded.  Moreover  we  require  that 
-A_,  fL  ,  -^ti..  etc.,  and  all  their  first  derivatives  vanish  on  the  boundr-ry  of  G,  while 
these  functions  possess  derivatives  of  all  orders  within  G,   These  conditions  on  u  are 
sufficiently  broad  to  apply  the  results  obtained  in  the  remaining  sections  to  the 
functions  u  consid'^red  in  connection  with  the  derivation  of  the  asymptotic  expansion. 

We  shall  now  derive  discontinuity  conditions  which  u  must  satisfy  on  any  surface 
on  which  it  is  discontinuous.   Let  G  therefore  be  a  region  in  which  u  is  discontinuous 
on  some  surface  ?5(x,  .•••,x  )=0  which  divides  G  into  two  domains  G^  and  G  .  We  shall 
apsn-ne  thnt  i   can  be  expressed  in  the  form  t-\jf  (x,  ,  •  •  •  ,x  .).   Let   I  ..  be  that  part 
of  r  which  bounds  G,  ,  and  let  P  ^  "be   that  part  ^Ich  bounds  G^. 

V'e  may  regard  u  as  consistinr  of  two  functions  u^  and  u  ,  such  that  u,  satisfies 
the  integral  equs.tion  in  G^  and  is  0  in  G  ,  and  conversely  for  u  .   Within  G^,  u^   there- 
fore satisfies  L(u)  =  0,  and  likewise  for  Ug.   In  G^  +  G  ,  u  =  u^  +  u  satisfies  (5.6). 
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Wp  dp^T^  conriltlnns   on    [n]    =  ii-i-"'ip  wki^h  hold   on   '6  =   0.      we  shpll  use   the  fact,    even 
for  d is continuous  u,    that  y   =  J     +  / 

Consider  first  a   typical  term  of    ("5.6)   for  i  and  j   =  1,2,-**,   n-1.     By  inte^^ra- 
tion  "by  parts 

/    u(a^'^A)_dv  =/      u(a><^il).(cbc.)  -  /  (a^'^il).u  .dv. 

Because  -fl-  and  J\  .  are  0  on  '  ■>   we  may  i!?;nore  that  part  of  the  "surface"  integral. 
In  the  integrand  of  the  remaining  surface  integral  we  must  understand  that  x  .  has 
been  replaced  by  the  exjjression  obtained  in  solving  ^  =  0  for  z.. 

By  another  integration  "by  parts  apilied  to  the  volume  integrs-l  we  obtain 

/u(a^^^A).^dv  =  /l(a^'^IV)  (dxj  -  /a^%u,(dxj  +/a^-5Au,  dv. 

let  us  note  now  that  the  direction  cosines  of  the  normal  to  ^  =  0  are 


±Jt7 


=  ±\  i^  . 


i     ' 


We  shall  supTJosp  that  the  plus  sign  in  front  of  the  radical  indicates  tho  direction 
of  the  normal  points  outside  of  G-  and  info  G  .  Also  let  a.,  i  =  l,2,«'',n,  denote  the 
direction  angles  of  this  "positive"  direction.   Let  us  also  use  the  notation  dw  for 
an  element  of  surface  on  ^  =  0.  '^e  may  rewrite  the  preceding  result  thus: 

/u(a^-^il),,dv  =  /Ii(a^%)A«5.dw  -  /^>-^riu  A0,dw  +  /a^^Au  ,dv. 

TTow  dw  =  dx,  •••d:c      ,    sec  a     =  dx;^«.«dx     ^   r-r-  .      Since  ur.  =        iH^     we   may  rewrite 
1  n-1  n  1  n-1  \0  j  *-      n 

the  ■preceding  result: 

/uCa^'^A)     dv  =  -/u(a^^il).vi,,(dx   )  -/a^-^A  uAi^.dw  +/a^%u,  ,dv. 

G  ^'^  0  ^    -^        "  0  \ 

In  the  first  surface  integi-al  we  must  remember  that  x  is  replaced  by  its  value 
w(x,,-'',x  _•]).   However  the  derivative  with  respect  to  x.  applies  only  to  the  x.  pre- 
sent before  the  substitution. 

To  achieve  integrands  in  which  only  A  and  Jl  .  appear  -  the  purpose  will  be 

t 

apparent  later  -  we  make  a  further  trensforma.tion  of  this  surface  integral.  We  note 
that 

^  (a^JA)  =  (a^-'A)^  +  ('^''■^)t  ^i  • 
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We  substitute  for  (a  '^■fi)     in  the  first  surface  integral  on  the  rifht  and  inte- 
grate by  parts  v/ith  respect  to  x.    the  integrand  containing  r —  (a  "^A).   The 
first  term  which  results  from  the  integration  "by  parts  must    "be   taken  over 
the  "boundpry  of  0,  which  boundary  lies  on  the  boundary  |   of  G.   Since  the  inte- 
grand contains  -A.  ,  which  vanishes  on  P  ,  this  integral  vanishes.   There  results 

/u(a^-'-n-).jdv  =/^^^-^^   (uwj)(dx^)  +/(a^JA)^un/^M/j(<i^^)  - /a'«^^UjX0^dw 

+  /a^-^-a  u,  .dv. 
%  ij 

We  now  make  simple  changes  to  obtain  1 

(').7)  y^u(a^^-a)^^dv  =/s^*^-^  |S-^t''l)V  ""^ij  (  ^^t^"   -^/(^'''^\^Vj''V* 

-  /'a^''nuA0,dw  +  /a^-^ilu.  ,dT. 

We  now  have  the  useful  form  of  those  terms  of  ('5. 6)  in  which  i  and  J  imn 
from  1  to  n-1.   We  consider  next  those  terms  in  which  the  derivatives  with  res- 
pect to  T.    and  T    (=t)  appear? 

/(ila^""),  udv  =/(a^''-n-).uX|!5,dw  -  /  (a^''il).u,dvi 

(5.8)     ^  ^ 

/  (iVa^'').^udv  =/(a^"il)^u\0^dw  -^^'^Au^\0^dw  +  / a^^'ilu.^dv. 

We   note   that   either  order  of  performing   the   integration  by  parts  with  respect   to 

X.   and  X     gives    the  sa:ne  result.  ^ 

in  ^^ 

We  consider  next   that   term  in   (5.6)  which   involves  — r  ; 

dx'^ 
n 

/  (a'^^il)     udv  =  /(a^^il)  u\0w  -  /  (a^'^il)  u  dv; 
\  nn  J^  n  ^^^  n  n 

(5.^)  /  (-''''A)     udv  =/(a^^A)  u\0,dw  -/a^"i\u  \0,dw  +/a"^Au     dv. 

./  nn  ^^  n       t  </^  n     t  (,  nn 

'^y^P.  fir^t-ortler   terms   in    (5.6)  are   treated  as   follows: 
(5.10)  -  /  (II b'^),  udv  =  -  /Ab^'uX^  dw  +  /Ab\dv. 

?'inally,    the   cAu  term  in    (5)   is   left  as   is.      Th?it   is, 
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(5.11)        /  ciludv  =/  ciludv. 


All  of  the  steps  from  (5.7)  to  (5.11)  could  "be  performed  with  respect  to  Gg 

instesd  of  G^ .     The  only  difference  is  th&t  \  would  136  replaced  "by  -X.  hecause  the 

normal  points  out  of  G     into  G  .  We  now  apply  (5.6)  rememherinf  that  y=  y  +  y  . 

2      1  ,  ,^  ,^ 

We  note  hpfore  collecting-  the  results  in  equations  (5.7)  to  (5.11)  that  the  vol- 

ume-integrel  terms  remaining  on  the  ri/?ht  side  add  ut)  to  +J^L{u).     Since  L(u) 

^1 
holds  within  G  and  all  quantities  in  L(u)  approach  finite  values  on  the  houndary 

of  G  this  volime  Integral  adds  uo  to  zero.  Similarly  f or  y  rLL(u).  "^rther, 

3ince  the  results  over  G.  and  G  are  the  same  except  for  sipn  we  may  comhine  them. 
HowcTer  those  VEdues  of  u  and  its  deriyatives  which  olitain  on  0  =  0  when  approached 
from  G  may  differ  from  triose  which  ohtain  when  approached  from  6  .  We  shall  use 
the  notption  [u]  ,  [u.]  ,  etc.,  to  indicate  the  differences  u^-u-,  '^T4~''^pit  etc. 

The  result  of  comhlninr  the  terms  in  (5.'')  to  ("5.11)  with  the  corresponding 
terms  for  G^  will  of  course  he  0  'because  of  (5.6).  Also  since  we  are  con- 
cerned in  the  paper  as  a  whole  with  the  case  in  which  the  coefficients  of  L(u) 

are  independent  of  x   =t,  we  shall  use  this  fact.  We  note  that  the  ternis  which 

n 

appear  on  the  right  sides  of  equations  (5.7)  to  (5.11)  contain  either  -H.  or  il-  . 
as  a  factor.  We  argue  that  -^  and  -H-.  may  he  made  0  on  0  =  0  independently  of 
«ach  other  while  retaining  the  otherwise  arbitrary  nature  of  the  -H- -function.   In 
view  of  this  fact  the  coefficient  of  -0.  and  the  coefficient  of  il   must  each  be 
zero  in  order  that  the  surface  integral  be  0.  We  consider  first  the  terms  in  -^  ^' 
These  all  contain  lul  as  a  factor.  If  |u|  Is  not  zero  then  we  find  that  on  0  =  0 

r.    B^'^.^vA    +  2  E  a^'^^.  +  a"X  =  0.' 
i,,1=l    ^   ^  ^  i=l     ^ 

We   recall   that   <f)  =  t-^u.      Fence  6=1,   0.=  -\v.  .     We   see   therefore   that   these   terms 

t  11 

amount    to 

n  .  . 

(5.12)  Yl       a     ?5  J     =0  on  0  =  0. 

i,-'=l  '    '' 

In  other  words,    the  di<^contlnuity  surface  0  =  0  must  satisfy  the  condition    (5.12) 
if  u  is   discontinx^ons   there.     Eauation   (5.12)    is  known  as   the  characteristic  condi- 
tion    .  T^anation   (5.12)   is   not  strictly  a  partial  differential  equation  since   it 
need  hold  only  on  ^  =  0.    Howerer  we  may  oonrert  It  to  a  proper  differential  equa- 
tion in  '1/  T^y  using  0  =   t-^ir  =  0. 

iF]     For  a   full  discussion  of  the  characteristic  condition  see  Coiarant-Hilhert,    loc. 
cit.,   v. II,    Chap.   6,  article  1. 
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We  now  consider  those  terms  in  equations  (5.7)  to  (5.11)  which  containJV,. 
Since  these  terms  must  elso  add  no   to  zero  we  may  write 


n-1  ^_  ,,  n-1   ,^_  n-1  n=il 

ij=l 

n-1 


(5.13)  21   Ha^^  V.a^t  +  T.      a^-'WM'i^t*  t  ^^J^i^l^  -  t  a*J[u  ]^ 

^^  ^   i,j=l     ^  -^   i,J=l  ^  ^  J  ^  i,1=l     J 


i=l     ^       ^    k=l      ^ 

Ve  recall  the  fact  that  ^  =  t-V  ,  and  hence  that  0^  =  1,  0^^=  0,0^^=  -  V^,  pj^j  =  -  \f^y 
i,.    =   0.  ''/e  use  these  facts  to  rewrite  (5.13)  as 

(5.1M  -  r  [-]-%,-  E  a'-^[u.i0  -.  E  ^'-^ w^v  E'  -'-'L-.ii^i 
i,j=i      ^^'  i,.i=i       -^  i.j=i     ^   ^  i,j=i     -^ 

-2E  a'"[u,]0,  -  a-[uj^-t  [u]h\  =  0. 
i=l       ^  k=l 

The  second  and  fourth  terms  comhine.      We  may  use  equation   (5.12)   to  replace   the   third 
term  "by   fu.]  ("ZjT'&^^^.i/).    -  a^^ii-i.).     The  result  thus  far  is 

(5.15)    i     Ma'*^^ii-2E     s>'5[u]^     -E  Mt\  =0       on  0  =   0. 

This  result  may  "be  called  the  first  discontinuity  condition  (aside  from  the 

characteristic  condition).   It  may  be  converted  into  an  ordinary  differential  equation 

in  rul  as  follov/s.  We  introduce  the  concept  of  differentiation  in  the  direction  trans- 

15 
verse  to  0  =  0.   The  direction  itself   is  given  by 

dx.   n   .,      ^ 


(5-16)    x^=  jr  =  5Ia% 


> 

k=l    "■    """i 

where  $  =  E  ^  "^i  ^^  • 
i,.i=l      ^ 

The  derivative  of  a  function  u(x^,'«-,x  )  in  this  transverse  direction  is  then  friven 
by  the  directional  derivative 


M^y^^i=r    a^J  1^  0.  =  E 
;^s  4-,  Bx,  i  ,'-7  ,    ax.  i^j  ,*-7. 


a  "u. 


i^i  ^=^i   ^   it7=i       ^^i   ^   itj=i      ^  J 

We  may  therefore  rewrite  the  first  discontinuity  condition  (5.15),  after  a  change  in 
sign,  as 


15.  See  CoTirant-Hilbert,  loc.cit.,  p. 35!^,  Since  the  normal  to  0=0  has  the  direction 
numbers  0   i=l,2,...,n,  the  transverse  direction  is  perpendicular  to  the  normal 
and  therefore  lies  in  the  surface  0=0. 
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(5.17)  [u]  (  i:  a^''^^^  +  X:  b\)  +  2^  =0  on  0  =  0. 

This  ordinary  differential  eqia^tion  will  he   called  the  (first)  transport  equation  , 
As  s  traverses  a  path  which  has  at  ench  point  the  direction  of  the  transverse  direc- 
tional derivative  (5.16),  equation  (5.1")  gives  the  behavior  of  the  discontinuity 
M  along  this  path.   It  will  he  noted  that  (5.1?)  is  an  ordinary  linear  homogeneous 
dif-ferential  equation.  The  equations  (5.16)  are  called  the  equations  of  the  char- 
acteristic rays  and  they  fire  the  rays  of  ordinnry  optics. 

We  may  ahhreviate  the  writing  of  (5.17)  hy  introduction  of  the  operator  L(u)  = 
n    .  .      n   ,  -^n 

Y2    a  '"'u   +  IZ  ^  '^'  "^^^^  (5.17)  reads  as  follows  ^: 
i,.i=l    ^-^   k=l 

(5.18)  L(0)  [u]  +  2  -^  =  0. 

1\'e  shall  specialize  the  characteristic  equation,  the  discontinuity  condition,  and  th« 

transport  equation  to  the  three-dimensional  wave  equation,  for   the  wave  equation 

2^2^2    22 
^1  •*■  ^2  •"  ^3  '  "^  ^t  =  ° 

and  for  P  =  t-v|r,   the  characteristic  equation   (5.12)  "becomes  the  true  partial  differ- 
ential equation 

(5.19)  I!    VK?  =  n^. 
i=l 

The  first  discontinuity  condition   (5.15)  "becomes 

(5.20)  -[u]    (11  TV,.)  +  2(  -Z.  [ti/lVi-n^[uJ  )  =  0, 

i=l  "       i=l 

wherein  t  is  replaced  "by  \^,      If   we  use  vector  notp.tion  and  renie:n"ber  tha.t  in  (5.20) 
t  is  to  "be  replaced  'by  \Uj  then  we  may  rewrite  (5.20)  as 

—  1ft 

(5.21)  /ikvbi]  +  2  gr8d[u]    •  grad  \i/  =  0. 

The  equation    (5.18)   for  the  rays  may  also  "be  specialized  to  yield 


16.    "Riad  we  not  assumed   that  the  coefficients   of  L(u)  =  0  were   independent  of  x     we  would 
have  obtained  from  the  derivation  in  this  section  the  sll^tly  more  general  transport 
equation  n  ..  n       ,  £_       <  i  ^W 


i,i=l        ^J    k=l       ^    i,j=l  "^     ^  J 


later  theory  in  this  paper  does  not  carry  over  to  this  more  general  case. 

17.  See  Courant-Hilbert,  loc.cit.,  p.36l. 

18.  This  result  was  obtained  by  Bremmer,  loc.cit.,  p.U21  and  also  Lowell,  loc.cit. , p. 282. 
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(5.22) 


^^i  =  "di"=-^i     ^  i  =  1,2.3  J 


dt     2 
di-=  -^    • 

One  sees  from  (5.22)  that  it  is  "both  possi^ble  and  convenient  to  replace  the  parameter 

s  aloae:  the  rays  "by  the  time  t.  Also^  since  the  cliaracteristic  surface  is  given  "by 

t=vj/,  we  see  from  (5.22)  tliat  the  raj^s  are  orthogonal  to  the  family  of  characteristic 

surfaces  in  3-3p6.ce  obtained  by  regarding  t  as  a  parameter. 

All  of  the  derivations  of  this  section  apply  as  well  to  the  non-homogeneous 

equation  (2.2)  as   to  (2.1).   If  one  considers  a  point  (x,t)  in  space  where  no  sources 

are  present,  that  is,  where  the  non-homogeneous  term  is  zero  or  even  continuous,  no 

19 
additional  terms  appear  in  the  characteristic  equation  or  discontinuity  condition  . 


19.  To  take  into  account  the  discontinuity  which  the  source  term  h(x)f (t)  introduces 
on  the  discontinuity  s\irface  0=^=0,  some  slightly  more  involved  considerations 
are  necessary.  See  Kline,  M.,  loc.cit.,  on  the  treatment  of  the  source  term. 
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VI.     Derivation  of  the  Higher  Discontinuity  Conditions  and  Higher  Transport 
Ea-ujRtions 

The  startin.'r  point   of  oiir  derivation  of  the  first  discontiniaity  condition  and 
first  transport  eqxintlon  was 


;.6)     J^MiCi 


(5.6)      /  ^MiCl)dv  =   0. 
G 

In  this  eq-a;?tion  Xt  Is  ar"bitra.ry  except  for  the  fact  that  it  and  its  first  deriva- 
tives vanish  on  the  boundary  P  of  G-  and  ^   must  possess  derivatives  of  all  finite 
orders  in  G.   In  view  of  the  arbitrariness  of  -ft  we  may  replace  it  hy  -ft-  ,  and  we 
now  utilize  the  additional  requirement  already  mentioned  in  Section  Y   that  the  first 
derivatives  of  il  .  vanish  on  I  .  We  find  therefore  that 


..1)    J. 


(6.1)     /  uM(n.  )dv  =  0. 


those  terms  which  contain  derivatives  with  respect  to  x.  and  x,,  i,j  =  l,2,...,n-l, 
in  addition  to  a  derivative  with  respect  to  t.  We  shall  treat  a  tj-pical  term 

(6.2)    /u(a^'^A).,,dv  =  /u(a^«5A)   (dx„)  -  /  u.  (ila^'^) .  .dv. 
^G^       tij    v/0       ij   n   J^^   t       ij 

In  the  surface  integral  it  is  understood  that  x  or  t  is  replaced  "by  v)/(x,  ,...,x  -) 
hut  the  subscripts  1,  j  imply  differentiation  which  would  obtain  before  the  substi- 
tut  i  on .   iTow 

We  use  this  relation  in  (6.2)  and  obtain 

(6.3)  /-(»'^'^)„/-y^fej(-'''-f^)ite„)  -^»(«'^-^)n''j(^„'  -  /^'^-'■')iJ*- 


We  apply  intef^ration  by  jyarts  to  the  first  surface  integral  on  the  right  side,  the 

t  to  X,.  The  first  l 

t.' 

u(a^''lL)^(dx  )(dx^). 


Integration  being  with  rpspect  to  x,.  The  first  term  which  results,  namely, 

t.' 


/ 


Bdy  of  0 

boimdary 
obtain  from  (6.3)  the  fact  that 


vanishes  because  the  boimdary  of  0  is  on  I  and  XL  and  XL  vanish  on  l   .   Fence  we 
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(6.4)  J^ni^^A)^^^,,  -   -Aa^-'tD^  |a-(a.^)-  A(«'^^)u''j<<^n>  "/"t  <''''")l/- 


1 

^^e  shall  use  next  the  relation 


Hence  from  (6.^), 


(aj  .A)j  =  g|-  (a'^A)  -  (b*JA  j^vj. 


We  may  integrate  hy  parts  with  respect  to  x.    in  the  first  integral  on  the  ri|?:ht  side 
and  use  the  fact  that  r\.  is   0  on  the  Tioimdary  of  f).     Then 

(6.5)^^u(.^^^A),,^d.  -/^^^''^)'^,T-^   (^^n)  *^(^''^>t-i  %K)- j-(^''^^it-/^^n) 

-/u.(rVa^'^)      dv. 

To  ohtain  terms  in  the  surface  integrals  which  involve  only  SX.   and  time-deriva- 
tives of  Jl.  we  use  the  relation 

in  the  third  surface  integral  on  the  right  side  of  (6.5).  We  integrate  hy  parts  the 
resulting  siirface  integral  which  contains  the  complete  partial  derivative  with  res- 
pect to  X.,  use  the  fact  that  A.  =  0  on  the  boijndary  of  0,  and  obtain  for  i,j  =  1,2,..., 
n-1. 

(6.6)^u(a^-^A),,^dv  =yJ(a^Jxi)^(dxJ  ^/^^^''^\^1  '^,^'\^  ^ f^^^' ' ^\^^^n^ 

■.^u(a^JXV)^^^^^^(d.^)  -^u^«la^J).jdv. 

We  now  consider  those  terms  in  (6.1)  which  involve  derivatives  with  reBpect  to 

X.  and  X  .  i=1.2 n-1,  and  a  second  derivative  with  res-nect  to  x  or  t.  We  treat 

in  n. 

a  typical  term,  namely 
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/^(^''^)itt  =  P^''^\t^^n^  -  L  ^-''^ht  V- 


Since 

G^  "0  1  p  G^ 

We  integrate  the  first  surface  integral  "by  parts  and  use  the  fact  that  fl      is   0 

Xr 

on  the  hounds ry  of  0.     Then 

(6.7)  /u(a"A)j^^  -  -/a"^t|^  te„)  -  A"j^„«ite„)  -/(-"^'„V'- 

G^  pip  Gj^ 

Vi'e  consider  next  the  term  in  (6.1)  which  involves  three  differentiations  with 
respect  to  x  : 

(6.8)  A(a"il),„  =/CaXt(^n)  "  /vXt^^" 

We  consider  next  the  second-ordet  terms  in  (6.1),  treating  first  those  involving 
diff crentintion  with  respect  to  x,  ,  k=l,2, . . . ,n-l,  and  x  ; 

-^u(h^JI)^^dv  =  -^(h^il)^u(dx^)+^  (h^A)^u^dv. 

Using   (h^-O-)-    =  |— (h^A)  -  h^-TL,   w,  ,  we  ohtain 

—    /u(h'"n),  .dv  =  -/  |-(h^Xl)^(d3,   )  +/h^A.    n/.u(dx   )  +/    (h^il)  u.dv. 
^Gj_  ^*  ^^0  ^\  n      y^         t     5c         n       V^_^  k  t 

Integrntion  hy  parts  in  the  first  surface  integral  plus  the  fact  that  A  is  0  on  the 
"boundary  of  i   yields 

(6.9)  -  /^.(.'^JLJ.^dv  ,yJ,''A|^(d.^_)  .J^^n^  V<^n>*/(''''^)l.^d,. 

As  for  the  second-order  terai  in  (6.1)  involving  two  differentiations  with  res- 
pect to  X  , 

(6.10)  -  /  u(b*i\)^^dv  =  -  /(h*A).u(dx  )  +  /  (b*A).u^dv. 

y^^  tt      v/^      t    n    J^^  t  t 

The  first-order  term  in    (6.1)   ie   treated  es   follows: 
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(6 


.11)     /A.cudv  =    /acu(dx  )  -/ ilcu.dv. 


We  intend  to  add  the  results  of  steps  (6.6)  to  (6.11)  to  obtain  an  expression  for 
J   TiM(Xl  )dv  =0.   We  then  intend  to  do  likewise  for   /  uM(rL  )dv  =  0  and  add  this 

resTilt  to  the  first  one.   This  sum  will  contain  on  the  right  side  -  /u.M(n)dv. 

y(j  t 

Now  this  voliime   inte^ra]    is  precisely   (5.6)  except   thst  u.    here  replaces  u  there. 
We  shall   therefore  use   the  resxilts   of  section  V   (equations    (5.12)  and   (5.18)    in  parti- 
cular)  to  write  thrt 

(6.12)  -yu^M(Ja)dv  =  -y  [u^.]At     fl    a^\^j  ^dw     +yiL(L(0)[u^]+   2  -^  )  Xdw 


i.j=l 


-  /AL(u^)dv. 


The  Tolume  int-grsl  on  the  right  side  must  he  understood  as  the  integral  ever  G^  pl\is 
the  integral  over  G  .  However,  since  u  is  required  to  have  continuous  third-order 
derivatives  off  0=0,  it  follows  by  differentiating  (2.1)  thpt  u.  satisfies  L(u)  =  0 
in  G-  pnd  in  G  .   Moreover,  the  discontinuities  in  u  and  its  derivatives  a.re  required 
to  he  Jump  discontinuities.   Hence  the  volume  integral  in  question  must  he  zero. 

Ve   shall  now  stun  the  results  of  eqioations  (6.6)  to  (6.11)  and  add  to  these  the 
analogous  equations  for  Gg.  For  the  resulting  volume  integral  on  the  left  side  we 
get  zero  hy  virtue  of  (6.1).  For  the  resulting  volume  integral  on  the  right  side  we  use 
(6.12).  Also  we  replace  the  (dx  )  in  (6.6)  to  (6.11)  hy  \0  dw  in  the  case  of  G,  and  hy 
-X0.dw  in  the  case  of  G  ,   Instead  of  writing  down  the  entire  surface  integral  which 
the  summs,tion  yields,  let  us  note  thst  the  integrand  contains  terms  in  fl,  -Tl  ,  and 
A  .  .  .  We  shall  treat  these  separately.  Neglecting  the  \dw,  we  get  for  the  coefficient 

[u]  I  ll  a^^  vu^  Wj  0^  -  2a^*  y^^  ^^   +  a^^t  ) 

If  we  use  the  fact  that  0  =  t-\v,   we  may  say  that 

(6.13)    ^tt->  w.r  ^^^  0i 0j  . 

We  now  collect  the  coefficients  of  -H.. .  These  are,  neglecting  \dw, 
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»t 


k=l  i»j=l 

Replacing  derivatives  of  y/  "by  derivatives  of  0  =  t-y  and  expanding  derivatives  of 
[u]  ,  wherein  it  is  -understood  that  u  =  u(x,,x  -...x  ,,\);),  we  obtain 

i,J=l     ^         i.J=l     "^    ^  i,J=l       ^''  i.j=l       ^  i,j=l         ^ 

i=l     ^  *     i=l     ^  ^  k=l      ^     ^i,j=l    ^  •' 

on  n    .  . 

If  we  use  the  fact   that  E  a  '(^^0^  =  0^  we  may  replace  the  second  and  seventh 

terms  "by  -a  "|jal0^,  and  the  resnltin^  expression  agrees  precisely  with  (5.1^). 
t~  t  t 

Fence  we  may  apply  the  stens  lepdin^  from  (5.1^)  to  (5.18)  and  write  that 


(6.1M  -\->-   [uJE  a'-5(i^«5  -  i(?»[u]  -  2 


eM 


the  first  term  coming  from  the  vol^ime  integral  (6.12)  and  the  remainder  from  other 
terms  in  (6.6)  to  (6.11). 

There  remains  now  the  coefficient  of  i\  .  These  terms,  aside  from  the  factors 
Xdw  and  ^^^are  as  follows: 

n-1   .  .  ;^2r  1    n-1  ,  Bful  sFu.! 

i,j=l      1  j   k=l     He 

The  last  two  terms  come  from  the  volume  integral  in  (6.12)  and  the  rest  from  other 

terms  in  (6.6)  to  (6.11).   In  (6.15)  it  must  be  rememherftd  that  \y\   =[u(x^,Xp, .  .  .x  _i»\;rll 

We  may  now  determine  how  the  a.uantities  in  (6.13).  (6.1^)  and  (6.15)  f're  related. 
If  [u]  j^  0,  then  the  argument  in  section  V  hased  on  the  independence  of  fL  and  XL, 
showed  thfit  the  charactrristic  equs-tion  and  the  first  transport  eqtiation  hold.  Then 
the  coefficients  of -H-^.  and  Xl,  as  given  in  (6.13)  and  (6.1^)  are  0.   If  ful  =  0  on 

20.  In  section  V  we  argued  that  the  coefficients  of  -H-and  A^  in  the  surface  integral 
are  independently  zero  hecause  of  the  arbitrariness  in  H.   Then  if  [u]  ^   0,  the 
characteristic  equation  holds  and  so  may  te  used  here.  However  we  shall  show  in 
?  moment  thPt  even  if  [ul  =  0  but  [u^"]  j:  q   then  the  characteristic  equation  still 
holds. 
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)j  =  0,   then  certainly  the  coefficient  of  ^.  .    here  la  0  and,   since  the  first 
transr^ort  equation  is   0,  we  may  throw  out  this  part  of  the  coefficient  of  i^.. 
A^flln  erPTiln^  from  the   independence  of  ^  and  XL    we  see  flret  from   (6.1^) 
thpt  if    1J5  J    ^  0  on  {5  =  0,   the  charscteristic  equation  muet  hold.     That  is, 
discontinnitles   of  u.   must  lie  on  a   surface  ()  =  0  which  satisfies   the  char- 
acteristic equation  even  if  u  itself  is   continuous  on  <5  =  0.      Secondly,   the 
coefficient  of  ^l    niust  also  he  zero.     This   last  fact   therefore  yields   the 
second  transiDort  equation  which  reads 

''i'his  equation  must  hold  on  0  =  t-^  =  0.  The  operators  on  the  right  side  apply  to  [ul  = 
[u(x-,...,x  ,,vi/3«  I't  will  ^e  noted  that  these  operators  involve  only  spatial  differ- 
enti?tions   or  differentiations   inside   or  along  0  =  0.     Eqxaation    (6.16)   is  a  non-horao- 

reneous      linear,   first-order     ordinary  differential  equation  in   ju.l   which  holas  along 

?1 
the   rays  '   , 

To  solve  equation  (6.16)  we  note  that  the  non-homogeneous  terra  woiild  "be  obtained 
from  the  solution  of  the  first  transport  equation.  This  solution  would  "be  ohtained  as 
a  fimction  of  s.  However,  eqimtion  (5.16)  gives  the  x.   es   f^jncuions  of  paramters  a. 
and  s.   This  family  of  rays  must  be  on  the  surface  0=0  and  hence  the  family  depends 
upon  n-1  paramters  and  s.   3y  inverting  this  system,  that  is,  solving  for  the  n-1 
parameters  and  s,  and  substituting  these  solutions  in  [tijjWe  obtain  this  quantity  as 
a  function  of  the  n-variables  X-.,...,x  ^and  then  may  obtain  the  derivrtives  which 
appear  on  the  right  side  as  functions  of  the  i. .  We  may  then  replace  the  x.  in  the 
derivatives  by  their  values  in  terms  of  s  to  solve  the  equation  (6.16). 

We  can  now  see  how  to  derive  the  transport  equations  for  Lu^J  and  the  higher 
time  discontinuities.  We  start,  by  nna.logy  with  equation  (6.1),  with  the  equation 

(6.17)  /  uM(A^^)  dv  =  0. 

We  now  sup^»ose  that  fl    .    and  all  its  first  derivatives  vanish  on  P.     Again  we  regard 

the  Integral   over  G  as  a  sum  of  integrals  over  G.   and  G  .  Analogously  to    (6.2)  we  may 
state  that  r  r  fl 

^^u(a^^A)^.^^dv=^u(a^^n).^^(d.J  -  J ^^^^^^^ ^) ,^^^-    y 

(6.18)  =  ^u(aiJa)^^^(dx^)  -^-t(a'^^^)i^(d-n)-/-tt(-'''^-)i^<^- 


21.    Compare   Courant-Hllbert,    loc.cit.,  pp. ?6 2- 36 3  and  equations    (12)   on  p. 363  in  parti- 
calpr.      The  qmntity  J  there   is   the   right  hand  side  here.  3qtiation    (6.16)   above   is, 
apart'from  the  method  of  deriv^-^tion,   a  more  general  result   than  equation   (12)   there. 
We  have   the  form  of   the   transport  equation  for  the  case  where   the  characteristic 
surface   is   any  0=0.      Moreover,    the   relationship  of    (6.16)   or  the  Courar.t-Hilbert 
eq-oation    (12)   to  the   transport  equation  for   [u"]    itself  is   now  established. 
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We  now  compare  the  Inst  two  terms  here  with  the  two  terms  on  the  right  side  of  (6.2), 
The  dlfferencpfl  are  in  sl^n  end  in  thnt  u.  here  replaces  n  there.  We  may  treat  the 
Inst  two  terms  here  ps  in  steps  (6.2)  to  (6.15)  «nd  ohtain  the  various  coefficients 
of  ri  ,  A^,  and  A,^^  piven  "by  (6.1?)  to  (6.15)  after  changes  in  sign  and  substi- 
tntion  of  u  for  n. 

The  first  term  on  the  rlpht  side  of  (6.18)  is  analofrous  to  the  first  term  on 
the  right  side  of  (6.2)  except  that  C\.      replaces  A  .  Again  we  may  repeat  the 
steps  from  (6.2)  to  (6.15),  noting,  however^  only  those  contrihutions  to  (6.13), 
(C.Ik)   snd  (6.15)  which  come  from  the  surface  integral  of  (6.?,)  and  that  the  coeffi- 
cients of  A^^,  fL^,  and  Pi   will  now  he  coefficients  of  -O-^tt  -^  tt,  -^  t  respect ivtly. 
The  following  will  result  for  the  coefficients  of  the  various  derivatives  of  -A-. 

(6.19)  i\,„  ->    K     t_  «'^^»J 

(6.20)  n,,   ->  -    [uj   Y:_  n'H.^.  -  L.(0)M    -   2  -^^ 

1 » (3"  J- 

n  .  .  ^  ^W        ^     ii     ^    M 

(5.21)   ^  ->+[uj  z.  -''Mj-*-  '^(^)  W  -^  2-^r-%^_,^    iT^ 

■"  i  .  ■1=1  1  » J"-'-  ■'•       J 


1 ,  ,1=1  1 .  J  -L 

n-1  ^  a[u] 


i7J=i       ^^i^^j  ■  fel         ^^ 


n-l        ..   o    \niA        n-i     ,       o  u. 

(6..2)    a  ->    n    »<.i  _£  .  ^  ,1^  ^ .  ,r„^]  .  i(0)t^j-2  ^^ 


We  now  may  argue  as  we  did  in  connection  with  the  derivation  of   (6.16)  from 

('^.1?)   to    (6.]5).      If   [uj   ^  0     on  0  =  Oj  the  c^amcteristic  equfition  and  first  and 

second  transport  equations  hold.     Hence  the  coefficients  of  A        ,  -f^..,  and    -^. 

ttt    tt        t 

vanish^ and  since  ^  is  srhltrary  we  get  the  third  transport  equation: 

If  |u|  »=  0,  "but  |u.  I  ^  0  on  0  =  0,  then  as  we  saw  above,  the  characteristic  equation 
holds^ end  of  course  the  first  transport  equation  holds  automatically.   Moreover, 
arguing  from  the  independence  of  H.  and  n.,we  see  from  (6.2l)  that  the  second  trans- 
port equation  holds.  From  (6.22)  we  obtain  (6.2?).   If  [u]  =  0  and  [uj  =  0  then 
the  first  and  second  transport  equations  hold  automatically.   If  at  the  same  time 
L'U'^^]  ^  0^  then  (6.21)  shows  that  the  characteristic  equation  holds,  and  we  have  the 
fact  that  discontinuities  of  [u  J  which  are  not  discontinuities  of  u  or  u.  must 
nevertheless  lie  on  a  ^  =  0  which  satisfies  the  characteristic  equation.  Moreover, 
(6.22)  shows  that  (6.23)  holds  in  this  case  too. 
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It  is  apparent  that  we  can  o'btain  higher  transport  equations  from  our  'baBic 
integral  equation.  Systematization  of  the  process  would  make  the  procedure  mech- 
anical. 

We  may  now  relate  the  general  theory  of  discontinuities  to  the  problem  of 
obtaining  an  asymptotic  solution  of  our  partial  differential  equation.  The  theory 
applies  to  any  discontinuous  solution  and  hence  certainly  applies  to  the  pulse 
solution  whose  discontinuities  furnish  the  coefficients  of  the  expansion.  What 
distinguishes  this  soliition  from  the  others  is  the  Initial  conditions  which  woiild 
uniquely  determine  the  solutions  of  the  transport  equations.  We  consider  this 
point  in  the  next  section. 


-3^- 

VII.  Initial  Conditions 

If  the  transnort  and  hif;her  transport  equations  are  to  be  used  to  furnish  the 
coefficients  of  the  aajTnptotic  expansionj Initial  conditions  are  needed.  No  general 
prescription  can  be  given.  Fowever^  it  dops  seem  possible  in  several  cases  to  state 

what  the  initial  conditions  should  be 

itac 

a)  If  a  plane  wave  of  the  form  e    impinges  on  a  perfectly  conductin/r  surface, 

e.g.,  the  exterior  of  a  parabolic  cylinder,  and  one  seeks  the  asymptotic  expansion 
of  the  reflected  field,  then  v;e  may  use  the  fact  that  the  sum  of  the  incident  and 
reflected  field  is  zero  on  the  conductor  and  hence  require  that  each  coefficient  in 

the  ES;yTnptotic  expansion  of  the  reflected  field  be  the  negative  of  that  of  the  inci- 

ikx 
dent  field.  Since  e    is  the  first  and  only  term  of  its  own  asymptotic  expansion, 

ikx 
in  whlcb  e    is  the  phase  factor  and  the  amplitude  coefficient  is  1,  we  may  require 

that  the  initial  value  of  the  first  term  of  the  reflected  field  be  -1  and  the  initial 

value  of  each  of  the  other  coefficients  be  0, 

b)  Supi-ose  the  theory  applicable  to  a  scalar  problem,  let  us  say,  with  a  point 

source  in  a  non-homogeneous  medium.  The  scalar  potential  of  a  Hertzian  dipole  is 

ikr 
in  a  homogeneous  medium.  This  term  is  also  its  own  asymrptotic  expa.nsion  in 

powers  of  X.  To  obtain  the  corresponding  potential  in  asymptotic  form  for  a  dipole 

source  in  a  non-homogeneous  medium  we  require  that  each  coefficient  of  the  desired  ex- 

ikr 
pansion  approach  at  r=0,  the  corresponding  coefficient  of  e .  To  avoid  the 

r 

infinities  Involved  we  note  first  that  the  analogue  of  r  in  non-homogeneous  media  is 

the  V^,  where  K  is  the  expansion  coefficient  of  the  reciprocal  of  curvature.  Hence 

we  multiply  the  n*-th  term  of  the  desired  expansion  by  (yZ)  a,nd  require  that  this 

i^r  n 

term  ermroach  the  n-th  term  of  the  expansion  of  e    multipled  by  r  .  Of  course  the 

r 
Iptter  term  will  be  zero  except  for  n=l. 


c)  Initial  conditions  may  be  obtained  from  the  "boundary  conditions  imposed  in 
problem  I  defined  in  section  II.  The  pulse  solution  U(x,t)  is  required  by  definition 
to  satisfy  the  boiindary  condition  (2.1c)  where  f(t)  is  Y\  (t).   Then  U(x,0+)=h(x)  and 
tj   (2.1a),  U(x,0-)  =  0.  Hence  [u(x,o)|]=  h(x).  Moreover,  for  x  on  s  ,  fu  (x  0)1  ••  Q 
and  likewise  for  the  higher  time  derivatives.   If  the  variable  t  is  used  as  the  para- 
meter along  the  rays,  then  we  do  have  the  Initial  values  for  the  transport  equations, 
at  least  insofar  as  the  field  which  propagates  along   rays  issuing  from  S.  is  concerned. 
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